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ABSTRACT 

We report here on a survey of N-body simulations of encounters between spherical galaxies. 
Initial systems are isotropic Jaffe models. Different sets of mass ratios, impact parameters 
and orbital energies are studied. Both merger remnants and systems perturbed after a non- 
merging encounter are analysed and compared to real-life elliptical galaxies. The properties 
of merger remnants show a large variety. Merger remnants resulting from head-on encounters 
are mainly non-rotating prolate spheroids. Merger remnants from models with J or b + are 
triaxial or mildly oblate spheroids supported in part by rotation. The velocity distributions are 
biased toward the radial direction in the prolate case and the tangential direction in the oblate 
case. Non-mergers are affected in various ways, depending on the orbital characteristics. We 
conclude that many of the global properties of real-life ellipticals can in principle be attributed 
to a merger of spherical progenitors. 

Key words: galaxies:interactions- kinematics and dynamics- structure- elliptical - numeri- 
cal simulation 



1 INTRODUCTION 



Elliptical galaxies, for a long time believed to be simple rotating 
spheroidal systems, are now known to be far more complex. From 
the study by Bertola & Capaccioli (1975) we know that support 
against gravity is provided by random motions as well as rota- 
tion. Remarkably the surface brightness distribution of ellipticals 
closely follows the r 1 ^ 4 law proposed by de Vaucouleurs (1958), 
although a more general r 1 ^" profile, with 1 < n < 10, (Sersic 1968) 
appears to provide a somewhat better fit (Caon et al. 1993, Bertin 
et al. 2002). Several dynamical models have been developed try- 
ing to reproduce this r 1 " shape (e.g. van Albada 1982). The shape 
of elliptical galaxies is in general triaxial and goes from oblate to 
prolate (de Zeeuw & Franx 1991). With few exceptions elliptical 
galaxies contain only small amounts of gas, and can be considered 
as dissipationless systems. 

At present there exist two main views on how elliptical galax- 
ies are formed. On the one hand a huge gas cloud may contract 
rapidly and form stars, resulting in an elliptical like galaxy. This is 
known as the gravitational collapse theory. On the other hand el- 
liptical galaxies might form via assembly of small building blocks: 
encounters of small systems and the subsequent mergers will give 
rise to larger galaxies. Elliptical galaxies formed through this pro- 
cess could even be the result of mergers of disc galaxies. This is 
known as the hierarchical merging theory. 

Since the early work by Toomre & Toomre (1972) pointing 
towards merging of spiral galaxies as a likely key factor in galaxy 
evolution, a large effort has been devoted not only to observe in- 
teractions and mergers of galaxies but also to model these pro- 



cesses and the resulting systems. On the observational side, features 
found in such detailed studies include: tidal tails, bridges, shells and 
counter-rotation (see Schweitzer 1986 and Schweitzer 1998 for re- 
views of the observations). Modelling of such features in individual 
systems with N-body simulations has in general been quite success- 
ful (see Barnes 1999 and Naab & Burkert 2001 for recent reviews, 
and for the various features: Hemquist and Quinn 1988, Merrifield 
and Kuijken 1998, Balcells and Quinn 1990, Hernquist and Barnes 
1991, Balcells and Gonzalez 1998, Toomre 1978, Higdon 1995, 
Weinberg 1998, Garcfa-Ruiz et al. 2001). 

Early papers on interactions of spheroidal systems, (van Al- 
bada & van Gorkom 1977, White 1978, White 1979, Farouki et 
al. 1983, Aguilar & White 1985) focus on the exchange of mass 
and energy. The main aim of these and of several later studies was 
to explore the conditions under which a merger will take place and 
to probe the general morphology of the merger remnants. For ex- 
ample, Navarro, in a series of papers (Navarro (1989), Navarro & 
Mosconi (1989) and Navarro (1990)), find that merger remnants 
resulting from an encounter of non-rotating de Vaucouleurs mod- 
els have a tail in the light distribution departing from the r 1 / 4 law. 
These simulations, with a few exceptions, have been done with 
one thousand particles at most. They therefore lack the resolution 
needed for a detailed study of the kinematics of the endproducts. 
In recent years mergers of spheroidal galaxies have not received 
much attention, although a number of studies continue along the 
lines of earlier work (Okamura et al. 1991, Vergne & Muzzio 1995, 
Combes et al. 1995, Levine & Aguilar 1996, Seguin & Dupraz 1996 
and Makino and Hut 1997). Detailed study of the structure of the 
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merger remnants, in particular the kinematics, is still largely unex- 
plored. In this respect the structure of merger remnants resulting 
from encounters between spheroidal galaxies is still poorly known. 

Until recently the merger rate of galaxies was crudely known. 
Therefore it has been difficult to put the results of the work de- 
scribed above into the context of galaxy evolution. As shown re- 
cently by van Dokkum et al. (1999) in their study of clusters 
at medium and high redshift, a large fraction of elliptical (or 
spheroidal) systems are undergoing mergers (up to 50 per cent). 
Therefore, the study of the end products of encounters between 
spheroidal galaxies is important for our understanding of ellipti- 
cal galaxies. This has led us to the following questions: (1) what is 
the relation between the input parameters describing an encounter 
(masses, impact parameter, relative velocity) and the parameters 
describing the merger remnant (shape, kinematics)? (2) how do the 
parameters of the merger remnants compare with those of real-life 
ellipticals? 

To answer these questions we have carried out simulations of 
collisions between spherical systems with a variety of input pa- 
rameters. We chose three input parameters to keep our experiments 
simple but at the same time reasonably complete. The models were 
let to evolve with different orbital conditions and with different 
mass ratios between the progenitor systems. The results are pre- 
sented in the form of a review of the main characteristics of the 
end products. Results from these simulations concerning the evolu- 
tion of the Fundamental Plane of Elliptical galaxies due to merging 
can be found in Gonzalez-Garcfa & van Albada (2003). In Paper 
II we include the effects of the presence of a dark matter halo and 
compare with the results obtained here. 

The end products of these merger simulations show a large 
variety of properties. Their characteristics can be traced to specific 
initial conditions. For some features the impact parameter is the 
key to the final configuration, for other features it is the initial mass 
ratio of the progenitor systems. Either way we find that one may 
attribute several characteristics of elliptical galaxies to collisions 
between 'simple', i.e. spherical progenitor galaxies. 

We conclude that the various characteristics of elliptical galax- 
ies as observed now are not detroyed by mergers. In fact mergers 
between spheroidal systems (as well as discs) may well be resposi- 
ble for many of these characteristics. 



2 MODELS 

Below we describe the setup of our experiments and the integration 
method. 



2.1 Initial conditions 

As initial conditions we use purely stellar non-rotating isotropic 
spherical systems described by Jaffe's law (Jaffe 1983). The imple- 
mentation used here was developed by Smulders & Balcells (un- 
published). The distribution function (DF) is derived from the po- 
tential and the density, and the particle positions and velocities are 
obtained from the DF. 

For a Jaffe model the potential is given by: 



P(r)- 
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<P(r) 



n \r + rj) 



(1) 



where G is Newton's constant of gravity, M is the total mass of 
the system and rj is the half mass radius. The corresponding mass 
density is: 
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and the mass inside radius r is: 
T -M. 



r + rj 

Further, 
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For a spherical system with an isotropic velocity distribution 
the DF is a function of energy only (Binney & Tremaine 1987). For 
Jaffe's model the potential-density pair has an analytical solution 
for the DF given by: 



/(£) = 



M 



[F_(V2^)- V2F_(V^)- yl2F+(V£) + F + (y[2S)], 



(5) 



where S = (-Erj/GM) and F ± (x) is Dawson's integral: F ± (x) = 
e f dxe ±x (Jaffe 1983). Particles are placed in a sphere with 
positions and velocities provided by the above DF. 

When constructing the models a cut-off radius is imposed to 
avoid having particles at exceedingly large radius. The cut-off ra- 
dius is equal to lOxrj, corresponding to the 91% mass radius of 
the Jaffe model. This results in a small change of the structural pa- 
rameters. By cutting off the outter 9% of the mass distribution the 
half-mass radius of the model is reduced to the 82% of the theoret- 
ical rj. 



2.2 Units 

Non-dimensional units are used throughout. We adopt G = 1 for 
Newton's constant of gravity. The theoretical half mass radius of 
the Jaffe model rj and the total mass M of the smaller system are 
also set equal to 1. Then, for a Jaffe model (v 2 > =1/2 and the half- 
mass crossing time T CI = 2rp/(v 2 } 1 / 2 = 2 V2. The models may be 
compared with real galaxies using, for example, the following scal- 
ing: 

[M] = Mj =4x 10 11 M , (6) 
[L] = n = 10 kpc, (7) 
[r] = 2.4xl0 7 yr. (8) 
With these, the velocity unit is: 

[v]=414km/s. (9) 

2.3 Method 

We have used Hernquist's (1987, 1990) version of the TREECODE 
on an Ultra-Sparc station where a typical run (with 2 x 10 4 par- 
ticles) takes of the order of 10 5 seconds for about 10 4 timesteps. 
Softening was set to 1/8 of the half mass radius of the smallest 
galaxy (R]/2 = 0.82 and e = 0.1). The tolerance parameter was set 
to = 0.8. We have done a test run with 8 = 0.6. No major differ- 
ences were noticed, while the computational time is increased with 
a factor 2. Quadrupole terms were included in the force calculation 
and the time step was set to 1/100 of the half mass crossing time. 
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Figure 1. A test run was made to test the stability of our initial models. In 
the top panel we give the evolution of the mass inside different radii, the 
top line gives the 99 per cent mass radius, while the bottom one gives the 
5 per cent mass radius. The lower panel gives the ratio between the half 
mass radius (R1/2) and the effective radius (i? c ). For a Jaffe model this ratio 
should be 0.74. 



2.4 Stability of initial models 

We have checked the stability of our input initial model for 28 time 
units, which correspond to 10 crossing times (see section l2~2l . We 
use a model with 10240 particles and a softening, e = 0.1. The test 
shows (see FigureQtop panel) that the system relaxes for about 4 
time units and remains stable thereafter. This initial relaxation is 
due to the presence of the particle softening in the code. 

The ratio of the effective radius (R e , measured as the radius 
enclosing half of the mass in projection) and the half mass radius 
{R\ll) remains close to the value 0.74 given by Dehnen (1993) for 
a pure Jaffe model; see FigureQbottom. 



2.5 Initial parameters 

We vary the following three input parameters: the initial orbital en- 
ergy, the impact parameter of the orbit and the mass ratio of the two 
galaxies. 

Let M, R and a = (v 2 ) 1 / 2 denote mass, radius and velocity 
dispersion. According to the virial theorem the velocity dispersion 



is given by tr 2 = aGM/R where a is a structure constant appro- 
priate for a Jaffe-law system. Models with different masses were 
constructed following the scaling relation between mass and radius 
given by Fish (1964), i.e.: 

Mi IR\ = M 2 IR\ = constant. (10) 
Homologous systems then also follow the scaling relation: 

a\lR l =cr 2 2 /R 2 . (11) 

To calculate the initial orbital energy we proceed as follows: 
the total energy of the system is given by the sum of the internal 
energy of each system and the relative orbital energy: 



En + E i2 +E olh 

1,1, 
= --M x o-\--U 1 o- L 1 - 



(12) 
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where En and E{ 2 are the internal energies of galaxies 1 and 2 and 
E OTD is the orbital energy. Consider the case where at infinity the 
total energy of the system is zero. Then, 
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2 * 1 2Mi+M 2 

where Voo(0) denotes the relative velocity at infinity for a system 
with zero total energy. 

Given M\, M 2 , a\ and cr 2 we can calculate the corresponding 
value of Vco(0). Using equations llOlll ll and ll3l the velocity at infin- 
ity can be expressed as a function of the velocity dispersion of the 
smallest system at the beginning, cr\ : 



vim : 



1 + * 

M 2 



1 + 

M 2 



-3/2\ 



(14) 



In order to represent a range of energies, we have chosen three 
different values for Voo, given by: V 2 = 0, 0.5V 2 (0), and V 2 (0) 
respectively. 

With this choice for the orbital energy, £" or b > 0, we have 
parabolic (Vco = 0) as well as hyperbolic encounters. Part of the 
orbital energy, both in parabolic and hyperbolic encounters, may 
be lost in the subsequent interactions of the two systems, leading 
to a bound system with E OI \, < 0, in those cases where the interac- 
tion is strong enough. So cases with E orD < 0, although not included 
explicitly as starting points, will appear as a consequence of inter- 
action. 

The impact parameter (D) of the orbit was chosen as follows: 
we are interested in head-on collisions, D = 0, as well as encounters 
with D t 0, but with D < R out , where ^ out is the radius including 
99% of the mass of the most massive galaxy. Impact parameters 
D > Rout would lead to a weak interaction which is hard to follow 
numerically. We selected three values of D. D = (head-on col- 
lision), D = Rout/2 and D = R ut- To define better the border line 
between merger and non-mergers, some models with mass ratio 
2: 1 and D + have also been run with initial values of D slighlty 
smaller than the scheme presented above. 

The parameters of the initial models are summarized in Table 
Q where the units used are model units. The first column gives 
the model name. The number indicates the mass ratio, the small 
letter denotes the impact parameter, with h for head-on (D = 0), o 
for off set (D = R out /2) and g for grazing encounters (D = R ou t)- 
The final capital letter denotes the energy of the orbit with P for 
parabolic, H for hyperbolic and Z for zero energy at infinity. The 
models with mass ratio 2:1 and slightly different D are labelled 
with an r. For some specific cases we have chosen slightly different 
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Table 1. Input parameters. The columns identify the run and give the mass ratio, impact parameter and 
orbital energy. 



Run 


Mi : M\ 


ImpactPar. 


^orb 


Run 


M2 : M[ 


ImpactPar. 


^orb 


IhP 


1 : 1 
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5:1 








IhH 


1 : 1 
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1 : 1 
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7.07 
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2gP 


2:1 


14.14 
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7:1 


26.46 
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2:1 


10 













IgHr 


2:1 


10 


0.239 










3hP 


3:1 








WhP 


10:1 








3hH 


3:1 





0.387 


\0hH 


10:1 





2.039 


3hZ 


3:1 





0.775 


WhZ 


10:1 





4.078 


3oP 


3:1 


8.66 





lOoP 


10:1 
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3oH 
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3gP 


3:1 


17.32 





\0gP 


10:1 


31.62 





3gH 


3:1 


17.32 


0.387 











initial conditions, these are given by a letter w in the name (runs 
IhW, IwP and 5wW) . 

For most runs we use 10240 particles per galaxy. To avoid 
heating due to large differences in particle masses we used 25600 
particles for the most massive galaxies (M2 > 7). Model IgP was 
repeated with 10 times more particles giving basically the same re- 
sults as can be seen in the following sections. This is model IgPG. 

Models leading to a merger were evolved for at least 8 to 10 
dynamical crossing times of the merged system after merging, to 
allow the system to relax (reach virialization). Conservation of en- 
ergy is good in all the runs; variations are less than 0.5 per cent. 



3 RESULTS 

The present sample consist of 49 runs; the results are summarized 
in Tableland are discussed below. The merger remnants do resem- 
ble real-life elliptical galaxies. Figure|2|shows their surface density 
profiles. These can be fitted by a de Vaucouleurs (1958) profile for 
about 10 magnitudes. The innermost points fall below a de Vau- 
couleurs profile, this result is suspect however because these points 
lie inside one softening radius (e 1 ^ 4 ~ 0.56). 

In the presentation below we focus on properties that can be 
compared with observations. 




0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 

R ./4 R t/4 

Figure 2. The surface brightness of the merger remnants versus R ' . In 
this plot a de Vaucouleurs (1958) profile is given by a straight line. Note the 
slight upturn with respect to a de Vaucouleurs profile in the outer regions. 
The dotted line indicates the softening lenght. 
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Figure 3. Evolution of systems in run \hP. This is a head-on collision between two equal mass galaxies. Numbers at the top of each frame show the time in 
computational units. The first encounter occurs around time 85. The two systems develop a faint 'halo' after this first encounter. Many particles in these haloes 
will eventually escape. The end result is a prolate cigar-shaped E5 system. 



3.1 Phenomenology 

Our models show a wide variety of features, not only in the final 
state but also during the collision process. In Figures[3]|4]and[5]we 
give some examples of the evolution seen in the various runs. 

The first two, Figures[3]and[4] provide examples of the evolu- 
tion of initial states leading to merger. Figure[3]gives the evolution 
of a model with equal mass components in a rectilinear head-on 
collision (run \hP). A large fraction of the orbital energy is lost in 
the first encounter and rapid merging follows. 

Next, Figure|4| shows the evolution of a system with mass ra- 
tio 2: 1 . The two galaxies are placed in a parabolic orbit with an im- 
pact parameter slightly smaller than the outer radius of the larger 
of the two galaxies (run 2gPr). After the first encounter at / = 60 
the smallest galaxy features a tidal tail ('plume'); it also produces 
a small bridge of particles. The merger is complete around t « 210. 

Finally, Figure [5] shows an example of non-merging systems. 



The masses of the initial galaxies are 5:1 (run 5wW). The systems 
are placed on a slightly non head-on hyperbolic orbit. This is not a 
typical case for non-merging systems in our sample, and we give it 
here as an illustration. The initial parameters are such that we are 
close to the border between mergers and non-mergers (see figure 
[6j. The most massive system is almost undamaged during the in- 
teraction. However, the small system loses particles that end up in 
the potential well of the bigger galaxy. Also a fairly narrow bridge 
between the two systems is formed. A small cloud of particles leads 
the small galaxy as it runs away from its bigger companion (at the 
lower right-hand side of the frames). 

The space of the input parameters is given as a 3-D plot in 
Figure [6] Here the three coordinate axes are £ or b, D and M2/M1. 
Circles indicate mergers and triangles non-mergers. All the systems 
studied here involve intense galaxy interactions, leading to several 
distinct features, also for those runs where merging does not occur. 

All runs (but one, see below) with £ or b = in Table|2|result in 
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Figure 4. Evolution of the systems in run 2gPr. This is an encounter between two systems with mass ratio 2:1, where the galaxies meet in a parabolic orbit 
with impact parameter slightly smaller than the outer radius of the largest galaxy (initially at bottom). The final system develops a rotating bar. 



a merger. In these simulations the systems lose part of the orbital 
energy in the first encounters, which results in 'heating' of the two 
galaxies. Subsequent passes through pericenter ultimately lead to 
merging. In all cases particles escape, carrying away some energy. 
We have considered a large impact parameter in only one case with 
£orb = 0- F° r this run (IwP) the impact parameter is larger than 
the sum of the radii of the two systems; ultimately it may result 
in a merger but the required time would be larger than a Hubble 
time. The same applies to all systems denoted with '??' in Table|2| 
Note that model 2gP takes a long time to merge (close to a Hubble 
time). This is why this model although considered as a merger is 
also marked with '??' in Table|2| The boundary between mergers 
and non-mergers in the plane E ox ^ = in Figure[6|has not been fully 
explored. 

Also some runs with energy E m ^ > lead to mergers. Those 
models have D = 0, that is, they are head-on collisions. For these 



the dynamical evolution is similar to that of the head-on collisions 
with £ or b = 0, but this is only true up to mass ratios 5:1. 

In summary, mergers from the present sample without a dark 
halo are confined to a limited region in the 3-D space defined by 
orbital energy, impact parameter and mass ratio (see Figure^. For 
hyperbolic energies we have mergers only for head-on collisions 
and mass ratios close to one. 



3.2 Morphology of the systems 

The various morphological features of the end products are sum- 
marized in Table[3] In the column 'Class.' we give a 'Hubble type' 
morphological classification for the merger remnants. To do so we 
have fitted the isodensity contours by ellipses and we have taken the 
mean projected ellipticity inside R e of the system from 100 random 
points of view (see Figure[8j. Table|3|gives the median value (note 
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Figure 5. Evolution in run 5w. The initial systems are placed on a grazing hyperbolic orbit. This simulation does not end in a merger. Note the plume of 
particles ahead and behind the smaller system. 



that in Table[3]a value for the isophotal ellipticity at R e for a given 
point of view is listed). The table also gives the kind of interaction. 
We considered the system to be a merger if more than half of the 
particles of each system are bound together in a single object as a 
result of the interaction. When this condition is not met we classify 
the interaction as follows. If more than 10 per cent of the particles 
of one system end up in the potential well of the other we call it 
'strong' interaction. When there is an exchange of less than 10 per 
cent of the particles, it is denoted as 'weak' interaction. 

The more flattened systems result from simulations with ap- 
proximately equal mass components, roughly parabolic orbit and 
impact parameter smaller than R out (see Figure Q as an example). 
There is a clear trend from E5 to El as we go to higher mass ratios, 
as expected. Frequently the end product has type E3 (see Figure[8j. 
In Figure [8] we give the distribution of Hubble types for the final 
systems as calculated from the one hundred points of view. The left 
panel shows results for merger models, while the right panel shows 



non-mergers. In both cases we find a large range in ellipticities, 
from E0 to E5 for merger models and even higher for non-mergers. 

For mergers resulting from systems with £orb > the least 
massive galaxy is usually highly distorted, some of the particles are 
trapped in the main potential well and shells are formed similar to 
those found by Hernquist & Quinn (1988). Figure|5|gives the radial 
velocity versus the radius for the particles in the smaller galaxy 
for run 3hH. Some of the particles from this system are placed on 
shells, visible as thin curved features. Although these shells are less 
obvious than those found by Hemquist & Quinn (1988) it should 
be noted here that these authors explored minor mergers between 
spheroids (mass ratios larger than 10:1) to produce shells. From 
figure[9|we find that shells are also produced in major mergers (3:1 
in this case). 

When looking at the non-merging systems we notice that the 
interaction is stronger when the masses are similar. However it 
seems that interaction features, such as bridges and tails, are more 
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Table 2. Final systems. Column (1) identifies the run and (2) gives the time where the run was stopped. (3) 
indicates whether the final result is a merger (Y) or not (N). Models which might lead to a merger but will take 
too long to follow computationally are denoted with '??'. (4) provides the l.o.s. ellipticity for the isophote at 
R e (e), and (5) and (6) give the axis ratios b/a, and c/a. In (7) the ratio between the rotational velocity and the 
central velocity dispersion is given. For those models not leading to merger two measurements are included, 
the first one is for the less massive system and the second for the more massive one. 
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prominent if the masses are not equal, e.g. MilM\ of order 3. For 
non-mergers with high orbital energy the most massive galaxy is al- 
most undisturbed after the interaction, contrary to the least massive 
one. Part of its particles end up in the inner parts of its neighbour's 
potential well, while some are expelled and form a plume or 'halo' 



around the escaping core (see Figure|5}- In models with D close to 
zero a bridge may form. 
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Table 3. Morphological classification of the end products. The classification 
is indicative only; it refers to the value of the median when measuring the 
type from 100 random points of view. 



Figure 6. The space of input parameters: orbital energy, impact parameter 
(D) and mass ratio (M2/M1). All models of our sample are shown; those 
leading to merger are depicted as filled circles, non-mergers are depicted as 
open triangles. 
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Figure 7. Contour plots of the projected density distribution for the remnant 
in run IhP. This system, resulting from a head-on collision, has a prolate 
cigar-shape structure. Top view along the major axis, bottom view along 
short axis. 
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3.3 Prolate and oblate systems 

We have measured the axial ratios b/a and c/a for each system, 
where a, b and c are the principal axes of the ellipsoid. To do so, we 
have calculated the eigenvalues from the moment of inertia tensor: 



hj = J P(r) [Sijr 2 - XjXj\ d 3 x. 



(15) 



For a homogeneous ellipsoid with axes 2a, 2b and 2c, the 
eigenvalues of the inertia tensor are: 
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Figure 8. Distribution of mean ellipticity inside R e , calculated from one 
hundred random points of view. Left panel: results for merger end-products. 
Merger remnants cover the entire observed range of ellipticities, from E0 to 
E5. Right panel: results for non-merging systems. Note that the tails of both 
distributions are affected by 'classification noise' . 
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Figure 10. Axis-ratios of end-products in 'tri-axiality space' (de Zeeuw 
& Franx 1991). Runs are identified as in Table [T1 Triangles depict non- 
mergers. 



prolate one (de Zeeuw and Franx 1991). Non-mergers are depicted 
as triangles and are nearly spherical as expected. 

For our merger remnants, we find that head-on collisions (D = 
0) give nearly prolate end products, while for runs with D 4- the 
end products are triaxial or oblate spheroids. Merger remnants re- 
sulting from systems with mass ratios different from one are closer 
to spherical. 

Our models achieve a flattening of up to 30 per cent of the 
initial system. It should be noted that the flattening achieved is not 
large, to some extend this may be caused by using the particles 
inside a sphere for the calculation of the axis ratios. Some elliptical 
galaxies are more flattened than the flattest systems produced in our 
runs. 



Figure 9. In run 3hH with mass ratio 3:1 a number of faint shells are 
formed. This figure shows the radial velocity versus the radius for the par- 
ticles belonging to the least massive galaxy, i.e. the progenitor of mass 1. 
The particles from this system are placed, via phase wrapping, in shells that 
can be seen in this figure as curved, relatively sharp features. 



Ei = (b 2 + c 2 )/5, E 2 = (a 2 + c 2 )/5, E 3 = (a 2 +b 2 )/5. (16) 

If we define the axes such that E\ < E 2 < E 3 , it follows that 
the axis ratios are: 



a 


jE 3 +E 2 


-Ei 




V £1 + £3 


-E 2 


a 


lE 3 +E 2 


-Ei 




\j Ei +E 2 


-E 3 



(17) 



(18) 



Axis ratios for the merger remnants in our simulations have 
been calculated with these equations using particles whose distance 
to the center of mass is less than the half mass radius. Results are 
shown in Figure [TOl The top right corner (1,1) represents a sphere, 
the lower right corner an oblate system and the lower left corner a 



3.4 Rotation 

We have placed a slit along the major axis of the system as seen 
in projection for 100 randomly chosen points of view and we have 
measured the streaming velocity as a fuction of distance to the cen- 
tre as well as the central velocity dispersion, cr . Thus we measure 
the streaming velocity as a radial velocity along the line of sight, 
in a way that is similar to observations of ellipticals. cr was mea- 
sured inside a radius of 0.2 x R e , i.e. about two softening radii (see 
also Gonzalez-Garcfa and van Albada 2003). In Figure [HUpp we 
plot the ratio of the maximal streaming velocity, V max , and cr Q as a 
function of l.o.s. ellipticity at^ e - Observational data are also given. 
Note that our merger remnants when seen from random directions 
are mainly small ellipticity systems with slow to mild rotation. 

The last column in Table [2] gives the ratio between the max- 
imum streaming velocity and <x from a particular point of view 
perpendicular to the angular momentum vector. (This point of view 
corresponds in general to maximal flattening, that is, one sees the 
intrinsic flattening of the system.) This quantity has been plotted 
versus e measured from the same point of view, in Figure \H\ bot- 
tom. In this plot the merger remnants cover a large part of the space 
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Figure 11. Streaming versus ellipticity for simulations compared with ob- 
servational data. Open circles are high-luminosity ellipticals, filled circles 
are low-luminosity ellipticals while crosses are bulges (from Davies et 
al. 1983). Top panel shows the cloud of points (small dots) for the merger 
remnants and non-mergers measured from one hundred points of view. Bot- 
tom panel shows the models as seen from a point of view along the y-axis 
of the initial configuration. From this point of view one recovers in gen- 
eral the maximal intrinsic flattening. Small symbols give the observational 
data, large symbols are the merger remnants and small triangles are the non- 
mergers. For details see text. 



occupied by the observational data. Non-mergers lie closer to the 
origin. Note that all equal mass mergers produce remnants with 
high intrinsic flattening. High values of V max /cr is this figure are 
associated with large impact parameters. 

3.5 Boxiness-disciness 

For the end products of merger simulations we have calculated the 
deviation, S(<f>), of the isodensity contours from pure ellipses and 
write: 

6((p) = S + ^a„cos(n</>) + ^^b„sm(n(p). (19) 

If the ellipse is a good description for the isophotes and it has 
been correctly fitted, S, a\ ,«2> a 3 and the b n should be small. Then, 
if the isophote is discy, 04 will be positive and if it is boxy, a\ will 
be negative (see Binney & Merrifield 1998). To increase the sig- 
nal to noise ratio for our results an average over 60 configurations 
was made, always calculating the isophotes for a projection paral- 
lel to the intermediate-axis. These 60 configurations are different 
'snapshots' obtained by evolving the final systems a bit further. 
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Figure 12. Mean radial variation of 04 after equilibrium has been reached 
(based on 60 snapshots along the intermediate axis) for runs identified in 
the upper left corner of each panel. Panels on the left show a boxy structure 
inside R e , those on the right a discy structure. Where two models are shown 
on the same frame the first one is given by filled circles and the second by 
open ones. For details see text. 



Figure 1121 shows the variation of a\ with radius for seven 
merger models. Left panels show models with mass ratios 1:1, 
parabolic orbit, and different impact parameters. For all those sys- 
tems the ellipsoid is boxy inside R e . These three systems have a 
variety of characteristics. The top panel refers to a non-rotating ra- 
dially anisotropic system, while the other two are anisotropic rota- 
tors (see section l3~6l . The bottom left panel also shows the results 
for the run lgPG with an average for 6 snapshots, the results are 
quite similar to run IgP. 

The right panels in Fieure fTH show the a\ parameter for four 
models with non-equal masses. They have in common that the im- 
pact parameter of their orbits is equal to half the radius of the larger 
initial system. In these cases the end products have a mildly discy 
structure. This disciness is more evident for the models with mass 
ratio 3: 1 (run 3oP, middle panel) and 5: 1 (run 5oP, bottom panel). 

We also find, in agreement with Stiavelli et al. (1991) and 
Governato et al. (1993) that boxiness-disciness is subject to 
changes in the appearance of the same system due to projection 
effects: the same system can be seen as discy or boxy depending 
on the point of view. Our 1 : 1 models are mainly developing boxy 
isophotes. 

It is important to notice that this process of collisionless non- 
equal mass merging with an impact parameter D j= gives on 
average discy isophotes (see Figure 1121 right panels). Although 
perhaps unexpected, it can be easily explained in terms of trans- 
fer of orbital angular momentum. The small system brings part 
of the orbital angular momentum to the inner regions of the rem- 
nants, transferring this into internal spin angular momentum. The 
smaller galaxy thus disrupts and preferentially populates a plane 
perpendicular to the angular momentum vector, giving rise to discy 
isophotes. In an ideal case one could think of this process as a way 
to form highly flattened spheroids inside dark matter haloes simply 
by the encounter of a spherical luminous system with a dark matter 
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Figure 13. Pattern velocity as a function of radius for the end product of 
run loP. A bar is present that extends beyond the half-mass radius, and has 
a sharp cut-off at twice the half-mass radius. 
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Figure 14. Streaming and figure rotation in the end product of run ioP. 
Triangles: circular velocity from the potential, open squares: streaming ve- 
locity, and diamonds: pattern speed. 



halo which in that case must dominate the mass distribution. (The 
merged system would then have a high M/L value). We find the 
same process acting in model 5gP, described in section 3.7. 



3.6 Figure rotation and velocity anisotropy 

With the algorithm outlined in section 3.4, we also calculated the 
eigenvectors of the inertia tensor. We used these to study the evo- 
lution of figure rotation for the merger models up to 20 dynamical 
times after merging. This has been done for the rotating systems 
only. As illustrated in Figure [131 the results show that, in addition 
to particle streaming there is a bar-like rotating figure extending to 
about two half-mass radii. The bar is more prominent for low mass 
ratio encounters (1:1, 2:1). 

To confirm this point we have plotted the various relevant 
properties versus radius in Figure fT4l The upper curve represents 
the circular velocity calculated from the potential. The lower two 
represent the streaming velocity from two points of view perpendic- 
ular to the angular momentum vector. The pattern speed is indicated 
by the diamonds. The corresponding angular velocity is nearly con- 
stant inside R/R\/2 = 2. 

The initial systems have isotropic velocity distributions. The 
end products of the merger process are often anisotropic however. 
A measure of this quantity is given by the anisotropy parameter, ji, 
defined as: 



where cr t is the tangential velocity dispersion: 

2 2 2 

<r, = a a + a, h . 



(20) 



(21) 



We have measured this anisotropy parameter for models with 
mass ratio 1:1 leading to a merger. The results are shown in fig. l 1 51 

Systems resulting from a head-on collision develop a radial 
anisotropy (j3 > 0). Models with non-zero impact parameter show 
two types of behaviour. For small impact parameters, therefore 
small orbital angular momentum, we find no net anisotropy. For 




Figure 15. The anisotropy parameter /? for the mergers with equal mass 
progenitors. The solid line gives the initial model, dotted lines are models 
with D = 0, dashed lines and dashed-dotted lines represent models with 
D # (loP and IgP respectively). The thick dashed-dotted line represents 
model IgPG. 



large impact parameter, i.e. large orbital angular momentum, the 
systems develop a tangential anisotropy, (J3 < 0). 

Thus, the shape of the orbit, that is, the initial orbital angular 
momentum, is reflected in the final orbital structure of the stars. 



3.7 Noteworthy features 

Two interesting cases are runs 5gP and 10gP. The systems in these 
runs will probably merge in the end, but the evolution time is so 
long that the simulation was stopped long before. At the first peri- 
centre passage a small fraction (around 2 per cent) of the least mas- 
sive galaxy was stripped off and these particles settled into a disc 
inside the more massive galaxy, forming a rotating disc (see Figure 
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Figure 16. In due time, simulation 5gP probably will end up in a merger. 
The figure shows an intermediate state when an interesting situation occurs; 
only particles from the least massive system are shown. Part of the particles 
have left their parent system. They now lie in the potential well of the more 
massive galaxy (not shown here, centered at (-10,0)) and form a rotating 
disc as seen in the top left inset. 



4 DISCUSSION 

We have presented results of merger simulations of spherical, one- 
component systems, varying three parameters: orbital energy, im- 
pact parameter and mass-ratio. Dark matter has been omitted at this 
stage. In a following paper (Paper II) we address the differences be- 
tween systems with and without dark matter haloes. 

Whether or not a merger takes place is determined by the ini- 
tial conditions and our 'merger space' is in principle similar to that 
given by Binney & Tremaine (1987, Figure 7-9). The structure of 
the merger remnants, which show a wide variety of morphological 
and kinematical characteristics, also reflects the values of the ini- 
tial parameters. We find prolate non-rotating ellipsoids, as well as 
oblate systems with considerable rotational support. Our principal 
conclusions are as follows: 

• The properties of the merger remnants match many of the 
global features of elliptical galaxies. 

• The initial orbital angular momentum defines the shape and 
the orbital structure of the remnant to a high degree. Head-on en- 
counters produce prolate spheroids; off-axis encounters tri-axial or 
oblate systems. 

• Tumbling systems or oblate systems containing a rotating bar 
may form in collisions with high impact parameter and mass-ratio 
close to 1. 

• Remnants resulting from nearly equal mass mergers can be 
appreciably flattened and are either non-rotating (consistent with 
luminous ellipticals) or flattened oblate rotators. 

• Off-axis encounters of ellipticals with mass ratio different 
from 1 can result in the formation of a rotating disc of stars in the 
inner parts of the main merger remnant (due to material stripped 
away from the smaller galaxy during the encounter). 

• Boxy as well as discy deviations in the isophotes are obtained. 
Discy deviations are found especially in models with impact pa- 
rameter D t 0. Disciness-boxiness is point of view dependent. 



• Shells may form in collisions involving systems with a large 
difference in mass, in accordance with Quinn (1984), Hernquist & 
Quinn (1988) and Hernquist & Spergel (1992). We also observe the 
formation of shells for mergers of mass ratio 3:1. 



Obviously, further study is needed of the probability of find- 
ing initial configurations as studied here in 'real-life' systems, for 
instance by studying configurations in cosmological simulations. 
This would give a measure of the likelihood of 'gas-free' merging 
for the formation of ellipticals as seen today. 

In fact, observational evidence pointing towards merging of 
spheroidal systems as an important formation mechanism for el- 
lipticals is mounting. Bell et al. (2003) study over 5000 early-type 
galaxies at redshifts up to 1 . They find that in their sample there are 
few bright blue galaxies that may fade into the luminous red galax- 
ies (early types) seen today, and argue that luminous red galax- 
ies (E's) could be the outcome of mergers between less luminous 
bulge-dominated galaxies. Van Dokkum & Ellis (2003) observe 
traces of star-formation for a sample of galaxies in the Northern 
Hubble Deep Field. The mass associated with these regions of star 
formation seems to be small however and they suggest the possi- 
bility of mergers between bulge-dominated, gas-poor systems as a 
likely explanation. Further, Khochfar & Burkert (2003), from semi- 
analytical modelling, find that the fraction of mergers of bulge dom- 
inated galaxies (early types) increases with time, independently of 
environment. They also conclude that at least 50% of present day 
E's are relics of major mergers between spheroidal systems. 

Elliptical galaxies may be divided into two broad groups con- 
cerning their kinematics and orbital structure. One group con- 
tains the low-luminosity elliptical galaxies that are rotationally sup- 
ported and have discy deviations of their isophotes. The second 
group contains the high-luminosity elliptical galaxies that are ra- 
dially anisotropic, slow rotators with boxy isophotes. Clearly, be- 
cause the merger remnant models can be scaled to any size, our 
merger simulations can not explain such relations; they lack an es- 
sential element. 

The issue boxyness vs. disciness is complicated and we will 
not speculate on the possible interpretation of the observed corre- 
lations (Nieto et al. 1994, Bender et al. 1994, Kormendy & Ben- 
der 1996, Mihos & Hernquist 1996, Faber et al. 1997, Naab et 
al. 1999). Here we want to emphasize that mergers can produce 
both boxy and discy isophotes. The crucial parameter is the im- 
pact parameter: off-axis encounters produce mainly (oblate) merger 
remnants with discy isophotes while head-on collisions produce 
prolate merger remnants with boxy isophotes. But these results 
are complicated by the point of view dependence of boxyness and 
disciness. 

For the observables explored here, the formation of elliptical 
galaxies may need no disc intervention at all. This result would ex- 
tend the conventional picture of hierarchical merging where from 
assembly of complicated (disc) systems we would get to simple 
(E0-E1) systems, as well as more highly flattened systems. The 
present study shows that from nearly spherical systems one may 
also get more complicated flattened spheroidal systems. 

Mergers between spheroidal systems are likely to happen and 
as shown here they do not erase the characteristics of E's. For exam- 
ple, the intrinsic shape of the remnant can be directly related to the 
impact parameter of the collision. In other words, prolateness and 
oblateness are probably relics of the orbital parameters. Mergers 
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keep a memory of the dynamics of the encounter and some charac- 
teriscis of ellipticals can be built during the merging stage. Mergers 
involving disks and gas introduce a broader range of features, but 
E+E mergers may provide a way of producing many (although not 
all) observed characteristics of E's that should not be overlooked. 
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